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Abstract —By relying on a wireless backhaul link, relay stations 
enhance the performance of cellular networks at low infras¬ 
tructure cost, but at the same time, they can aggravate the 
interference issue. In this paper, we analyze for several relay 
coding schemes the maximum energy gain provided by a relay, 
taking into account the additional relay-generated interference 
to neighboring cells. First, we define spatial areas for relaying 
efficiency in log-normal shadowing environments and propose 
three easily-computahle and tractable models. These models allow 
the prediction of 1) the prohahility of energy-efficient relaying, 2) 
the spatial distrihution of energy consumption within a cell and 
3) the average Interference generated hy relays. Second, we define 
a new performance metric that jointly captures hoth aspects of 
energy and interference, and characterize the optimal number 
and location of relays. These results are obtainable with signif¬ 
icantly lower complexity and execution time when applying the 
proposed models as compared to system simulations. We highlight 
that energy-efficient relay deployment does not necessarily lead 
to interference reduction and conversely, an interference-aware 
deployment is suboptimal in the energy consumption. We then 
propose a map showing the optimal utilization of relay coding 
schemes across a cell. This map combines two-hop relaying and 
energy-optimized partial decode-forward as a function of their 
respective circuitry consumption. Such a combination not only 
alleviates the interference issue, but also leads to a reduction in 
the number of relays required for the same performance. 

Keywords - energy efficiency; interference; relay de¬ 
ployment; cellular network; decode-forward; relay coding 
schemes; log-normal shadowing; models for performance 
analysis 


I. Introduction 

In the urge to limit the energy consumption of cellular 
networks while maintaining service quality and ubiquitous 
access, relaying is a flexible and economical solution to 
enhance performance, eliminate coverage dead zones or al¬ 
leviate traffic hot zones. It is envisioned as part of next- 
generation cellular networks, along with pico- and femtocells 
jin . Unlike small cells, relay stations are not connected to 
the core network through a wireline backhaul connection 
but have to rely on wireless transmission to access the base 
station. This offers significant infrastructure cost reduction and 
deployment flexibility but, at the same time, can aggravate the 
interference issue. Exploring energy-optimized relaying jointly 
with interference reduction and choice of coding scheme opens 
new perspectives for efficient relay deployment. 

This work is supported in part by the FQRNT and the NSERC - 
Ultra Electronics Industrial Chair in Wireless Emergency and Tactical 
Communications. F. Parzysz and F. Gagnon are with Ecole de Tech- 
nologie Superieure, Montreal, Canada. M. Vu is with Tufts University, 
Medford, USA. (Fanny.Parzysz@lacime.etsmtl.ca; Mai.Vu@tufts.edu; Fran¬ 
cois.Gagnon@etsmtl.ca) 


A. Motivation and Prior work 


As highlighted in ||2l, the relay location within a cell can 
significantly affect the system performance. Its optimization 
is necessary to guarantee maximal gains within reasonable 
deployment cost and to avoid poor relay locations with neg¬ 
ligible gains. Substantial efforts have been paid in optimizing 
the relay location with regards to capacity coverage 

nil and energy |@, |^. The serving area of relay stations 
has generally been envisaged as a small circular zone around 
each relay ifl^ or covers the cell edge exclusively as in ||^. 
However, such serving areas are neither energy- nor capacity- 
optimized. An analysis of the serving area has been provided 
in d for capacity and in ifl^ for energy, via the newly- 
defined relay efficiency area. These works, however, account 
only for the transmit energy but not circuitry consumption, 
and consider only path-loss but not shadowing, which is a 
significant cause of signal degradation. 

On the one side, analyzing only the useful transmit power, 
i.e. the power radiated by the antenna, allows fair character¬ 
ization of the network upper-bound performance, as done in 
M. However, it is generally not sufficient for an energy- 
efficient analysis. Depending on the considered technology 
and hardware quality, the energy dissipated in circuitry for site 
cooling, network maintenance and signal processing (even if 
no transmission is operated), may dominate the overall con¬ 
sumption 00 such that the performance limits predicted 
by theory may not be realized. This is particularly true when 
the network does not operate at full load. 

On the other side, considering the transmit energy alone 
provides a basis for complementary analysis of multi-cell 
networks, as it is directly related to inter-cell interference 
(ICI), another great challenge of cellular networks. Unlike 
base stations, relay stations are generally equipped with omni¬ 
directional antennas and can drastically increase the interfer¬ 
ence to a neighboring cell. Several questions remain open: 
Does energy-efficient relay deployment necessarily lead to 
interference reduction? Does deploying few relay stations 
far from the cell edge but serving a large part of the cell 
generate less interference than deploying many relay stations 
transmitting at low power but potentially closer to cell edge? 
Answering these questions necessitates the understanding of 
whether there exists a trade-off between energy-efficient and 
interference-aware relay deployment in cellular networks. 

The ICI constraint has been investigated in the context of 
the relay placement for capacity and coverage enhancement 
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e.g. around a circle centered at their serving base station. 
In this case, interference appears to be very pessimistic, 
since each neighboring node (among the base stations and 
relays of each cell) is assumed to interfere with the reference 
user in each time slot with full transmit power. An exist¬ 
ing model for interference evaluation allows for interference 
from nearby stations to be computed in closed-form, while 
afar interferers are modeled as a continuum rather than a 
discrete set ||3l. While such models are particularly suitable 
for capacity maximization or coverage extension, they are not 
useful in analyzing or optimizing energy consumption because 
of the assumption that all stations transmit at the maximum 
power while the actual transmitted power can be much less. 
Understanding how energy consumption is distributed across 
the cell is essential for an energy approach and, to the best of 
our knowledge, no interference analysis has been performed 
for energy-optimized relay-assisted networks. 


B. Main contributions and Paper overview 


Our objective is to investigate how the deployment of 
decode-forward-based (DF) relay stations, specifically their 
number and location, can reduce the overall energy consump¬ 
tion of a cell and how it affects at the same time the perfor¬ 
mance of neighboring cells, due to the additional interference. 
First, we propose spatial definitions of relaying efficiency, 
namely the Efficiency Areas, as the cell area for which a 
given performance requirement is satisfied. For shadowing 
environments, we define: 


the Relay Efficiency Area (REA), inside which a user has 
at least a probability to be served by the relay; this 
REA extends the model in lll2ll to include shadowing, 
the Energy Efficiency Area (EEA), for which the energy 
consumption does not exceed E^; 


Eor the energy consumption, we account for (a) the transmit 
energy, (b) the additional energy consumed by the relay station 
for signal processing (decoding and re-encoding), (c) the 
energy loss at the RE amplifier and (d) the transmission- 
independent energy offset dissipated for network maintenance 
and site cooling, also called idle energy. 

Simulating the network performance offers wide possibil¬ 
ities but is time-consuming and can rapidly turns out to 
be unfeasible in shadowing environments. Moreover, they 
provide neither generalization to other simulation settings 
nor performance limits like an analysis based on capacity 
bounds. Therefore, we propose in this paper easily-computable 
models for relaying probability (REA) and energy consump¬ 
tion (EEA), which allow meaningful performance analysis 
without requiring extensive simulations. Such models have 
wide application and offer valuable support for load balancing, 
resource management or base-station switch off. 

Second, based on the proposed models for REA and EEA, 
we define a new framework for interference analysis. We 
characterize the average relay-generated interference to a given 
neighboring user and propose a model for it. Next, we define a 
new performance metric, called F, which balances the energy 
gain provided by relays and the additional interference they 
generate. 


Third, we apply the proposed models to investigate efficient 
relay deployment. With an energy-efficiency perspective, we 
characterize how the circuitry energy consumption, dissipated 
for signal processing and network maintenance, affects the 
network performance, as a function of the cell radius and the 
number of relays. We also highlight that energy-efficient relay 
deployment does not necessarily lead to interference reduction 
and conversely, an interference-aware relay deployment is 
suboptimal in terms of energy reduction. 

Finally, most approaches to efficient relaying are limited 
to the analysis of the simple two-hop relaying scheme. This 
scheme provides some performance enhancement but fails 
to capture the true potential of relays. Thus, in addition to 
two-hop relaying, we consider two energy-optimized partial 
decode-forward schemes of tH, which respectively minimize 
the network consumption and the relay consumption only. 
The second scheme is particularly relevant for analyzing 
relay-generated interference. Given the increased circuitry 
consumption of such optimized schemes, we draw a map 
showing the cell areas for the optimal use of each scheme that 
maximizes the energy-to-interference ratio F. The resulting 
map not only alleviates the interference issue, but also allows 
the reduction of the number of relays necessary to reach the 
same performance. 

This paper is organized as follows. The cell configuration 
and channel model are described in Section m The reference 
coding schemes and related model for energy consumption 
are described in Section Uni The models for REA and EEA 
are proposed and analyzed respectively in Sections ITVl and IV] 
The new framework for interference analysis is characterized 
in Section IVTl The application of these model in energy- and 
interference-aware relay deployment is covered in Section IVlIl 
and Section rVlIIl concludes this paper. 

II. System Model for a Relay-Aided Cell 

We present in this section the main assumptions on the 
cell environment and channel model. We use the following 
notation: BS stands for the base station, RS for the relay station 
and U for the user. Calligraphic TZ refers to the user rate. 
Upper-case letters denote constant distances, such as D, R for 
radius, or H for height. Lower-case letters stand for variable 
distances or angles, e.g. the mobile user coordinates {x, y) or 
(r, 6). Subscripts d, b and ^ respectively refer to the link from 
user to BS (direct link), from relay to BS and from user to RS. 
P is used for probabilities and E for expectation. Einally, DTx 
refers to direct transmission, while RTx stands for relay-aided 
transmission, in the considered coding scheme. 

A. Cell topology 

We consider an hexagonal cell with edge distance D;,, as 
depicted in Eigure [T] It consists of 3 sectorized base stations 
BSi, i S {1)3}, located above surrounding buildings. We 
assume a typical radiation pattern for each base station and 
consider the antenna gain as given in If^ . Each 120°-sector 
i is served by Nr relay stations, all equipped with omni¬ 
directional antennas. A given relay is at a distance of from 
its assigned BS. Einally, we assume a mobile user positioned 
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at {x,y). It is associated with the closest base station and the 
relay that provides the highest energy gain. We do not consider 
techniques for relay selection or multi-relays schemes but the 
proposed analysis can be adapted to such configurations. We 
define the maximal coverage by the outage requirement Pout 
that has to be maintained throughout the whole cell. 

B. Description of the relay channel 

We consider half-duplex relaying performed, without loss 
of generality, in time division. We assume that the relay 
operates on the same frequency resource as the user it serves 
(namely in-band relaying in LTE-systems). The multiple ac¬ 
cess strategy allows orthogonality between users within the 
cell, such that only one user is served for a given time and 
frequency resource. A downlink transmission is carried out in 
two phases of equal duration. We assume Gaussian channels 
with independent additive white Gaussian Noises (AWGN) 
with equal variance N on all links. We respectively denote hd, 
hi, and hr the channel coefficients from base station to user 
(direct link), from base station to relay (wireless backhaul link) 
and from relay to user. In addition, hi refers to the channel 
from the interfering relay to a user in a neighboring cell. 

We assume that the transmitted signal is degraded by both 
path-loss and shadowing. Note that we do not consider small- 
scale fading, since fading coefficients are rarely known for 
power allocation and terminals are designed to be sufficiently 
robust against such small-scale parameters, especially for 
relatively low mobility environments. Hence, this work is 
based on a per-block capacity analysis and can be understood 
as the mean relaying probability and energy consumption over 
a sufficient period of time such that small-scale fading is 
averaged out. 

The BS, RS and user have different heights, moving neigh¬ 
borhoods and transmission ranges, such that different links 
have different properties, particularly in terms of path-loss. 
To fit this heterogeneity, we consider the channel model 
proposed in the WINNER II project El, where the path- 
loss of link k G {d, 6 , r,/}, denoted 7 ^, is specified by four 


parameters Ak, B^, Ck and Dk depending on the global 
location of the transmitter and receiver (street level, rooftop...). 
The shadowing coefficient Sk is modelled as a log-normal 
random variable of variance of.. We assume that all Sfc’s are 
mutually independent. Channel gains are given by: 


with 


= 


Sk 

Ik 


Sk 

Kkd°‘>‘ 


{ oik = Ak/IQ 

10 logio i^k) = Bk + Ck logio (^) 

-\-Dk logio ((Htx - 1)(Hrx 


( 1 ) 

1 )) 


where Hxx and Hrx are the respective heights (in meters) of 
transmitter and receiver, d is the distance (in meters) between 
them and fc the carrier frequency (in GHz). 


HI. Coding schemes and Models 
FOR Energy Consumption 
We discuss the relaying schemes considered for analysis 
and their respective overall energy consumption. To simplify 
notation and focus on energy, we consider normalized trans¬ 
missions of unitary length, thus setting up a direct relation 
between energy and power. 


A. On the overall energy consumption 

The time-varying transmit energy, i.e. the energy radiated 
at the output of the RE antenna front-end, is referred as 
£:(RF) We assume individual energy constraints over the two 
transmission phases for the BS and RS, denoted as and 
respectively (with for maximum). Most literature 
only accounts for this useful transmit energy, which is a fair 
assumption for capacity or coverage analysis. However, it is 
generally not sufficient for an energy analysis. 

Eirst, significant energy is dissipated in circuit electronics 
for data transmission, especially in the RE amplifier which 
usually performs with considerably low efficienc y. W e con¬ 
sider the simplified, yet meaningful, approach of 1191], where 
the amplifier inefficiency is assumed linear in the transmit en¬ 
ergy and is characterized by a multiplicative coefficient, 

denoted rjn for the relay and rjE for the base station. 

Second, we account for the circuitry consumption related 
to a transmission, i.e. to signal processing at the encoder 
and decoder. As in 12], we model this consumption by the 
following energy offsets: Eg^’'^ for the BS transmission, 
and Eg^^^ for the RS reception and transmission, and for 
the reception at the mobile user. Such energy offsets mostly 
depend on the quality of electronics and on the complexity of 
the signal processing performed at terminals, notably given by 
the considered coding scheme. Using simple relaying schemes 
can help decrease this energy consumption, but at the cost of 
potentially degraded network performance. 

Third, in addition to the energy related to the transmission 
itself, we consider the transmission-independent consumption, 
also referred as the idle energy. This is the energy dissipated 
for site cooling, network maintenance and additional signaling, 
which is consumed whether or not data is transmitted. As 
proposed in El, it is modelled by an offset, consumed at 
each station and denoted as Eg'"^^ for relays, and Eg'^'^^ for 
base stations. 
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B. Coding schemes considered for analysis 

We now present the reference coding schemes and express 
for each scheme the overall energy consumption. In this work, 
we consider Gaussian signaling which accurately approxi¬ 
mates OFDM-based communications as used in practical 
systems including LTE and WiMAX. We consider downlink 
transmissions and assume that channel coefficients are known 
at transmitters to achieve the best possible performance. The 
power allocation aims at minimizing the energy consumption, 
given a hxed user rate TZ. 

1) Direct transmission (DTx): Data is transmitted directly 
from the BS to the user over the two transmission phases, 
and not just only during the hrst one as done in most 
literature. This two-phase transmission allows fair performance 
comparison since, in this case, both direct and relay-aided 
transmissions have the same delay and consume the same time 
resource. The energy consumption is given by: 

Edtx -f (2) 

Setting Nr = 0 gives the consumption of a reference scenario, 
where no relay station is deployed. However, when Nr > 0, 
the energy consumption of DTx should account for the idle 
energy dissipated at relay stations, even if DTx does 

not actually use those relays. 

2) Two-hop relaying (2Hop): Two-hop relaying is the sim¬ 
plest decode-forward scheme and thus, gives lower-bounds of 
the performance that can be achieved with DF-based relaying. 
In this scheme, the user sends its message to the relay station 
with rate 2TZ during Phase 1. Then, the relay decodes the 
message, re-encodes it with rate 272. and forwards it to the 
destination during Phase 2. The base station hnally decodes 
using only the signal received from the relay station (the direct 
link is ignored). The energy consumption is expressed as: 

E2H0P 

+ -P . (3) 


We denote E. 


(dsp) 


= E, 


(Rx) . p(Rx) 


-2Hop “ "-R + ^ additional energy dissi¬ 

pated at the RS for decoding and re-encoding, in comparison 
with the energy offset of DTx described in Eq. (|2|l. 

3) Optimized partial decode-forward schemes: Einally, we 
consider for simulations the partial decode-forward (PDE) 
schemes optimized for energy proposed in IB- In this, only 
part of the initial message is relayed, the rest being sent via 
the direct link. The base station divides the data message into 
two parts rrir and using rate splitting. During the hrst 
transmission phase, irir is broadcast to both the relay station 
and the user. At the end of this phase, only the relay decodes 
nir and then re-encodes it. During the second phase, the relay 
sends rdr and the base station jointly sends {mr,md) using 
superposition coding. At the end of phase 2, the user jointly 
decodes nir and rrid to recover the initial message. Several 
power allocations have been proposed in Olhll . In this work, 
we consider the two following schemes: 

• Energy-optimized PDE (EO-PDE), which minimizes the 
total transmit energy consumption and is both energy- 
and rate-optimal (referred to as G-EE in Olhll ). 


• Interference-at-Relay PDE (IR-PDE), which minimizes 
the energy transmitted by the relay only, thus minimizing 
the relay-generated interference (namely R-EE in lllhll '). 
We refer the reader to Olhll for the detailed power and rate 
allocation for nir and vrid- Such optimized schemes arguably 
require complex implementation and hne synchronization but 
they provide theoretical upper-bounds of the performance 
achievable with DE-based relaying. 

Little information is available on the circuitry energy con¬ 
sumed by such coding techniques. To account for the increased 
complexity of these schemes, we consider an additional energy 


offset E 


(dsp+) 
pdf 


to the overall consumption as follows: 


Eeo = -f -F prE^^^ 

+ + (e^"^) + NrE^^'^^) + E(4) 

In the performance analysis of Section rviii we will consider 
several values for . As for two-hop relaying, we denote 

-F + Epjf’’''’^ representing the additional 
energy dissipated in circuitry to perform this relaying scheme. 


IV. Relaying Probability and Relay Eeeiciency 
Area with shadowing 

We hrst characterize the relaying efficiency in shadowing 
environments and dehne the Relay Efficiency Area (REA) as 
the cell area where RTx has more probability to save energy 
compared to DTx. We analyze how to compute such area and 
propose a simplihed model for fast and accurate performance 
evaluation. Such analysis is described in detail since it will 
be used as basis for the energy and interference analysis, 
presented in Sections |V] and IVII 


A. A new definition of the Relay Efficiency Area 

Eor a given channel realization {hd, h},, hr), the area covered 
by a relay-aided cell can be divided into two geographical 
regions, depending on whether DTx or RTx should be per¬ 
formed. The wider is the cell area served by the relay, the more 
efficient can a relay station be considered. As highlighted in 
the introduction, various models are considered in the literature 
to characterize the serving area of a relay, but none are energy- 
optimized or account for shadowing. Here, we extend the 
dehnition of the pathloss-only Relay Efficiency Area (REA) 
in id to account also for shadowing environment. We dehne 
the REA as the set of all user locations for which relaying is 
statistically more energy-efficient than DTx or is necessary 
to satisfy a given outage requirement Pout- Both coverage 
extension and energy gains are included in the REA dehnition. 

Dehnition 1. The Relay Efficiency Area (REA) of a network 
in shadowing environment is dehned by the pair (^r,Pt). 
Any mobile user M (x, y) within the geographical area 
is served by the relay station with at least the probability Pt, 
either because RTx is more energy-efficient or because DTx is 
not feasible, i.e. Ar — {M {x,y) s.t. Pt < Prtx}, where 
Prtx {x,y) is the probability for RTx at user M {x,y). 

The REA depends on the user rate, the relaying coding 
scheme and the channel radio propagation. As an example. 
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Fig. 2; Relay Efficiency Area: simulation and model 
(Two-hop relaying, Nr = 1, D^. = 700m) 

the REA for two-hop relaying is illustrated in Eigure |2] for 
various and The dotted lines refer to the model 

for the REA proposed in Section IIV-BI Similarly, we can 
define Pdtx {x, y) as the probability for DTx at user M {x, y), 
either because DTx is more energy-efficient or because RTx 
is not feasible. Then, the outage condition is expressed as 
1 - Pout < Prtx {x, y) + Pdtx {x, y) for user M {x, y). 

B. Characterization of the relaying probability Prfx (a:, y) 

Eor readability, the relaying probability is analyzed for 
downlink two-hop relaying only, but is valid for other sce¬ 
narios, as discussed in Section IIV-EI We first focus on the 
case = 0. 

1) Notation for the energy consumption: We define Ed as 
the RE transmit energy consumed by the BS when DTx is 
used. Similarly, Eh and Er stand for the energy consumed by 
the BS and RS respectively when RTx is used. Considering 
Gaussian signaling, we deduce the energy consumption based 
on capacity formulas. To send data to user M {x, y), we have: 

Ed = ( 2 ^ - 1 ) 

Eh = (22^ - 1) iV2:^ (5) 

Er = (2^^- 1) Af^^, 

where N is the variance of the AWGN, Sk is the shadowing 
coefficient and jk the path-loss, as given in Eq. O- Then, we 
use superscript to refer to the non-shadowing case, when 
only path-loss is considered (i.e. the variance of the shadowing 
coefficient Sk is cr^ = 0). Assuming log-normal shadowing 
environment, we have 

Ek = log A/" {pk, crl) with pk = In (e^^A . (6) 

2) Probability for energy-efficient relaying: Eor a given 
channel realization, a transmission is relayed if DTx is not 
feasible or if RTx is more energy-efficient. This implies: 

Prtx {x, y) = PcR {x, y) + Per (x, y) (7) 


where Pcr and Per are as defined below. As we will see 
in next section, these probabilities are necessary to define a 
model for the energy consumption and interference. 

Definition 2. When the user-BS link is weak and data cannot 
be sent using DTx given the channel realization, relaying is 
performed to extend the cell coverage and maintain the outage 
requirement. This is referred to as the Coverage Condition for 
Relaying (CR), which occurs with probability Pcr. Given the 
energy constraints at BS and at RS, we have 

CcR={Ed>El^^ n Eh<El^^ n Er<E^^^}. 

Since energy consumption in each scenario is often inde¬ 
pendent, Pcr can generally be computed as 

Pcr = P[Ed> P (Sf, < P (s, < . 

Similarly, we define the Coverage Condition for Direct trans¬ 
mission (CD) for which data can be sent only via the direct 
link. It occurs with probability Pcd and is generally computed 
in closed-form as for Pcr. 

Definition 3. When both DTx and RTx are feasible, relaying 
is performed if it is more energy-efficient. This defines the 
Energy-Efficient Condition for Relaying (ER), expressed as: 

Cer = {Sk < E^'"^ n Er< E^"^^ n Eh + Er<Ed]. 

with k G {d,b}. The ER-condition has probability Per. 

Similarly, we define the Energy-Efficient Condition for 
Direct transmission (ED) for which DTx is more energy- 
efficient. It occurs with probability Ped: 

Ped = P [E^dM < n Er< - Per. (8) 

Thus, the outage condition at user M (x, y) is given by 
1 — Pout < Per + Pcr + Ped + Pcd 

< P (E{d,h} < n E;, < + Pcr + Pcd 

Due to power constraints. Per is obtainable by a triple 
integral over Sd, Sh and Sr, as in Eq. (|2^ of Appendix A. 
Thus, contrary to Pcr, it is not separable and may not exist 
in closed form. Even tough numerical computation can be 
envisaged using mathematical software, this approach rapidly 
becomes unsuitable for large networks, even considering the 
simple two-hop scheme. 

C. Proposed model for the REA 

To address the key issue of computation, we propose a 
closed-form lower bound for Per. 

Lemma 1. The probability Per for energy-efficient relaying 
is lower-bounded by the sum Plow = Plow + ^^^low where Pj^jJ 
and p[^2 ^6 given in Eq. (|9l) and (ITOl) at the top of 
page. Contrary to Per, the lower-bound Plow consisting of 
elementary probabilities that can be individually computed or 
have closed-form approximations. 

Proof: See Appendix lAl ■ 
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P™ = max (o, P (Ei+r < Ea) - [ P < Ei < + e[,""’) P i^Eb+r < + E^'"^) + P (e^'") + E^,'"^ < ) (9) 

pS =P <Ea< P (Et,+, < (10) 


Also note that, Pio* gives an upper-bound for Ped: 

Fed < P (E{aM < n - Pio*. (H) 

Model 1 (REA). The Relay Efficiency Area, as characterized 
in Definition [T] is modeled in log-normal outdoors propagation 
environments by the pair (.4 r,Pt), where 

.Ar = {M (x, y) s.t. Pt < Plow (a;, y) + Pcr {x, y)} . 

Plow and Pcr are given by Lemma [T] and Definition respec¬ 
tively. 

Any mobile user M {x,y) circulating within (yfRjPr) is 
also within (yfR,PT), such that it has at least the probability 
Pt to save energy via relaying. 

D. Accounting for the circuitry consumption 

Up to now, we have focused on the case = 0, i.e. ac¬ 

counting only for the transmit energy. These results, however, 
can be generalized to show the overall energy consumption. 
In this case, data is relayed if the overall energy consumed 
using RTx is less than the consumption using DTx. As defined 
in Section IIII-BI the overall consumption includes the RF 
amplifier efficiency (given by rjn and tjb) and the additional 
energy dissipated for decoding and re-encoding at the relay 
station (given by The energy-efficient condition for 

relaying becomes: 

= {k + e: + E(^=p) < E2 n n 

E° < I with E^ - log A/" {pk + In (pfc), al) 

Model 2 (REA, extended). The REA accounting for the 
overall energy consumption is modelled by (AIr,Pt), with 

-4r = {M(t,j/) s.t. Pt < Piow(a:,2/)+PcR(a:,y)}- 

Pj°^ and P^R are given by Lemma [T] and Definition |2] respec¬ 
tively, but with the following replacement: 

Ek —)> E^, ke{d,b,r} 

E,+r^ E^ + E° + E(^^P^ 

p(m) _^ p(m) (12) 

This model is illustrated in Eigure |2] for = 50mJ. 

E. Discussion 

1) Model validity: The proposed model has been validated 
under several outdoors environment settings and for both up¬ 
link and downlink transmissions. Eor the uplink, we generally 
have e[j'^^ < eL^^ for the user and relay energy constraints, 
implying that = 0. This model can be also used for any 


DF schemes for which the consumed energies Ed, E^ and 
Ed are log-normally distributed, such as the repetition-coded 
full DF scheme in where the user decodes data using 
maximum ratio combining on the signal received from the BS 
and RS during both phases. 

2) Model utilization: The primary use of the REA is 
to compute the spatial distribution of transmit energy con¬ 
sumption, as described in the next section. Thanks to the 
proposed model, many relay configurations and propagation 
environments can be analyzed in a reasonable time, which 
even allows finding the optimal relay location by exhaustive 
search. Moreover, this model can help decide between DTx 
and RTx when only statistics of the channel realizations are 
known at the transmitters (partial CSIT). In this case, the 
proposed model specifies for each user location the path which 
has the highest probability to save energy. 

Next, the concept of Efficiency Area has wide application 
since it is based on the network geometry and ensures a 
minimum performance. Such a framework allows the de¬ 
ployment of relay stations in an efficient manner, simply by 
locating relays such that hotspots or cell regions with poor 
performance are included within the corresponding Efficiency 
Area. This concept can be particularly applicable to non- 
uniform random user locations, e.g. a hotspot-type distribution. 
More specifically, the probability for energy-efficient relaying 
offers valuable support for network resource management, 
relay selection, load-balancing, scheduling or BS switch off. 
Eor example, switch-off can be decided if the probability 
to reduce energy consumption via relaying is above a given 
threshold Pt, i.e. if the considered cell area is included in 
(AIr, Pt)- 

V. Energy Consumption and Energy Efeiciency 
Area with shadowing 

We analyze the spatial distribution of the energy consump¬ 
tion within the cell and define the Energy Efficiency Area 
(EEA) as the cell area for which the average transmit energy 
consumption does not exceed Et. The EEA is related to 
the maximum energy necessary to transmit data at a given 
rate. Such a metric is relevant for performance analysis and 
encompasses fairness between the served users. Indeed, a relay 
station cannot be energy-efficient for users located in a certain 
area of the cell only, while showing unacceptably high energy 
consumption elsewhere. 

A. Definition of the EEA 

The average transmit energy, denoted as E con¬ 

sumed to send data to a given user M {x,y), is equal to the 
sum of the energy consumed by BS only when DTx is used. 
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Fig. 3: Energy Efficiency Area: simulation and model 
(Two-hop relaying, Nr = 1, D^. = 700m) 


and by both BS and RS when RTx is used. It is expressed as 
follows, with Eb+r = Eb + Er'. 


E 


£;(RF) 


= PcR E [Eb+r I Ccr] + IPer E[E;b+.| ^^er] 
+ PcD E [Ed I Ccd] + IPed E[E;d| ^^ed] ■ 


(13) 


Definition 4. The Energy Efficiency Area (EEA) of a network 
in shadowing environment defines the energy range across 
the cell. It is characterized by the pair [Ae^ Et): the average 
transmit energy consumed by any mobile user M (x, y) within 
Ae does not exceed Et, i.e. 


Ae = 



s.t. 


E 


Ei^n 



Eigure [3 illustrates Ae and its cotTesponding model for 
various values of E^. Also note that the areas for relaying 
probability and for energy Ae are distinct. Next, we 
describe the proposed model for the EEA. 


B. Analysis of the EEA with shadowing 

We focus on the computation of E [£'('^'^1] in Eq. (fTSl l 
for log-normal shadowing environments. Eor the coverage 
conditions CD and CR, both expectations are computed in 
closed-form: 

r E[£,+,|Ccr]= E 

< (14) 

[ E[Ed\CcD]=g(.d,E[^^) 


where the function g is given by: 


g{k, E) =exp[g,k + Y 


2 \ $ 


ln(E)-/jfc \ 

o-fc J 

<j) ^ ln(E)-/Xfc ^ 



(15) 


On the contrary, the conditional expectations E [Ed \ Ced] and 
E [Eg + Er I Cer] for energy efficiency can only be expressed 
in an integral form. We thus propose to bound both of them, as 
done for the probabilities Per and Ped in the previous section. 

Lemma 2. When RTx is more energy-efficient than DTx, 
the average consumed energy PerE [Eb+r \ Cer] in log-normal 
shadowing environments is lower-bounded by Ej™l, where 


p(ER) 

■^low 


= exp 


^Bb+r 


^b+r 

2 


Ip(l.ER) 

■^low 


+ g{b + r, 


p(m) 



fl ER') fl") 

Here, PjE ’ is computed similarly to P|E but considering the 
scaled distribution exp{(j‘lj^f)Eb+r rather than Eb+r- We recall 
that P{^2 and P{^2 are given by Eq. (|9]l and (fTOl i respectively. 

Proof: See Appendix IbI ■ 

Lemma 3. When DTx is more energy-efficient than RTx, the 
average energy consumption is equal to PedE [Ed \ Ced] and 
is upper-bounded by Eq. (fThl l at the top of page. is 

computed similarly to but with the scaled distribution 
exp(cr^)£'£; instead of Ed- 

Proof: The proof is similar to that of Lemma |2] ■ 

We now apply Lemmas |2] and |3 to deduce a model for the 
RE energy consumption. 

Model 3 (EEA). We model the Energy Efficiency Area (EEA) 
of a relay-aided network in shadowing environment by the pair 
{Ae, Et). The average transmit energy to any mobile user 
M (x, y) within .4 e can be assumed to be below Et, i.e. 


Ae = 


{x,y) s.t. E 


£(RF) 


< Et 


where E 


£(RF) 


and E(f). 


is given by Eq. (fTSl l but with the bounds 


Simulations show that E 


£(RF) 


tight upper-bound for E [iJ^RR^] and 


can be considered as a 
any user located with the 


model area Ae can be assumed to be within .4 e as well. The 
proposed model for the EEA can be extended to account for 
the circuitry energy consumption by using the same variable 
replacement as for the REA, given in Eq. (fTTI) . 




































VI. A NEW FRAMEWORK FOR JOINT ANALYSIS OF 
RELAY-GENERATED ICI AND ENERGY 
The assumption that all interferers are transmitting at maxi¬ 
mum power is not realistic within an energy-efficiency context, 
and we propose to characterize the impact of a RS on the 
interference imposed on a neighboring user, given the actual 
relay energy consumption. Such an analysis provides help¬ 
ful support for interference management techniques as well. 
Indeed, techniques to mitigate interference from neighboring 
cells, as considered for next generation cellular OFDMA sys¬ 
tems sia, require a fine understanding of the interference 
profile over time and frequency. This is particularly true in a 
network where both DTx and RTx occur simultaneously. 


A. Approximation of the relay-generated interference 

Based on the previous analysis for probability of energy- 
efficient relaying and overall energy consumption, we compute 
the average interference received at a mobile user located 
in another cell, at distance di from the interfering relay, 
as illustrated in Figure [T] To isolate the impact of a relay 
station, we do not consider other source of interference. This 
assumption is fair given that relays are usually equipped with 
omnidirectional antennas, as opposed to base stations. 

We denote Mj{xi,yi) as a given neighboring user, 7 / = 
Kid’^^ as the path-loss between the relay station and this 
user, and sj as the corresponding shadowing coefficient. The 
average interference received at Mi{xi,yi) is expressed as 


I{xi,yi) = Es 


E 




fl 

j jj 


= E 




exp {aj/2) 

Kidf 


(17) 


where E 


E.. 


(RF) 


is the energy radiated by the relay, averaged 


over Sd, Ss and s^. It is given by 


E 




= PcrE 




+ Per P 


I Cf 


(18) 


First, the energy radiated by the relay when DTx is not 
feasible, due to fading, is computed as: 

=PcR 5 (r',E^'")) 


PcrE 


(19) 


with g given in Eq. (fTsT l. Second, using a proof similar 

by the 


E 


(RF) 


C, 


ER 


to Lemma |2] the energy radiated PerE 
RS when RTx is more energy-efficient than DTx, is lower- 
bounded by with 

^ 2 ' 


Eiol = exp ( y 


11 . 

2 




( 20 ) 


where and are computed similarly to and Pj^^ 
in Eq. (|9]l and (fTOl l. but considering the scaled distribution 
exp(tT^)£^j. instead of Er- 


Model 4 (Approximation for I{xi, yi)). The relay-generated 
interference at a given neighboring user Mi{xi,yi), located 
at a distance dj from the relay station, is lower-bounded by 


I{xi,yi) = E 



exp (cr|/2) 
Kjd^^ 


= (PcRp(r,E("VEL) 


exp (cr|/2) 


where g is given by Eq. (fTST i and E®^ by Eq. ( l20l i. 

One can argue that upper-bounding the interference 
I{xi,yi) would be more suitable for performance analysis. 
However, we highlight that the proposed lower-bound is tight, 
as shown in Section IVII-AI It is thereby accurate enough to 
investigate interference-aware relay deployment. The proposed 
model for the ICI can be extended to account for the circuitry 
energy consumption by using the variable replacement of Eq. 
(fTSl i. similarly to REA and EEA. 


B. A new metric for analyzing energy and interference 

A relay station can provide significant energy gain and 
coverage extension for the cell it serves. But, at the same 
time, it is an additional source of interference, implying that 
neighboring cells experience an energy loss to maintain the 
same data rate for their own users. In consequence, a relay 
deployment is efficient if the achieved energy gain, referred as 
wcain, is higher than the resulted energy loss, referred as wloss- 
We propose to use their ratio as a metric to jointly capture the 
aspects of energy and interference. 

To evaluate the energy gain vcam, we consider a user 
M(x,y) served by BSi. We compare the energy 
consumed to send data to this user when BSi is not supported 
by relay stations (Nr = 0 ) and the energy consumed 

when BSi is supported by Nr relay stations. We have: 

< [a;, y] = [x, y] -f [x, y] 

[ -f -f -f -f -f . 

Here, the various energy offsets accounted in Eidie are de¬ 
scribed in Section IIII-BI We recall that we can focus only 
on Sector 1 since the three base stations BSi, BS 2 and BS 3 
within the considered cell are orthogonal and do use the same 
resource (in time and frequency). 

Similarly, to evaluate the energy loss bloss, we consider a 
user M{x,y), which is located in a neighboring cell i f I 
and performs DTx. We compare the energy consumed 

by the neighboring BSi to send data to this user when BSi 
is not supported by relay stations (Nr = 0 ) and the energy 
consumed to maintain the same rate, when BSi is 
supported by Nr relay stations generating interference. To 
isolate the impact of the Nr relays, we assume an ideal 
network, without any other source of interference. Denoting 
I (x, y) the interference received at M(x^ y) as defined in Eq. 
dnii of previous interference analysis, we have Mi 

{ [x,y] = geE^^^) [x,y] + -f 

I [a;, y] = [x, y] -f 

[ + j + 

Remark: In addition to energy gain, relay stations also 
provide coverage extension. To account for such extension and 
compute Bcain and bloss, we do not consider a power constraint 
for DTx. 
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Energy 

p(idle) 

pfdsp) 

'^2Hop 

riB 

p(Tx) 

_Eb_ 

IJ 

25mJ 

0-50mJ 

2.66 

p(Rx) 

+ '^11_ 

p(ni) 

(idle) 

■-R 

p(dsp+) 

■^pdf 

Vr 

90 

500mJ 

lOmJ 

0-50mJ 

3.1 

TlJ 



N 

-93dBm 


Pout 

0.02 


o-d 

6dB 


n 

3bit/ch.use 

Channel 

0-6 

3dB 

Others 

fc 

2.6GHz 


CT 7* 

4dB 


Hb 

30m 


o-i 

6dB 


Hr 

20m 


Normalized Tx (Is) 


Hu 

1.5m 


TABLE I: Simulation parameters 


Definition 5. To capture both the energy and interference 
aspects, we define the ratio T as 


t^Gain . , 

r = — with 

t^Loss 


t^Gain ^ 
'GLoss — 


{Nr=0)_p(Nr) 


r(Nr = 0) 


EEi 

i^l 


(Nr=0) 


EE; 

i^l 


(N,) 


y-' ^(Nr— 0 ) 


( 21 ) 


Here, t>Gain is averaged over all users served by BSi and 
wloss is averaged over all users located in the neighboring cells 
2 to 7, as depicted in Figure [T] If T > 1, the considered relay 
configuration is efficient, if 0 < T <1, the relay stations result 
in more energy loss for neighboring cells than they actually 
provide energy gain in their own cell. If T <0, relaying does 
not provide any energy gain, due to the circuitry consumption. 


VII. Performance analysis for Energy- and 
ICI-EFF iciENT Relay Deployment 

In this section, we first validate the proposed models for 
relaying probability, energy consumption and interference. 
Then, we jointly analyze the energy consumption and the 
generated ICI using two-hop relaying. In the last subsection, 
the impact of the relay coding scheme on the network perfor¬ 
mance is explored. If not specified, we consider the simulation 
parameters of Table H] taken from E El E m. For the 
channel gains, the direct link hd and interference link hj are 
modelled by scenario C2 of the WINNER II project lEsll . the 
RS-to-user link hr by scenario Bl, the BS-to-RS link hb by 
B5c. We recall that we consider normalized transmissions of 
unitary length, setting up a direct relation between energy and 
power. 


with = |m (a;, y) e Ar n M {x, y) ^ ^r| 

U |m {x,y) ^ Ar r\ M {x,y) G ^r| 

Similarly, the error ratio (jg for the EEA refers to the 
proportion of erroneous users, for which E < Ey and 


E 


e;(RF) 


> Et (or reversely), for some given threshold E^. 
For the interference analysis, we focus on the approximation 
of the average energy radiated by the relay and define the error 


ratio as the proportion of users for which E 


e: 


(RE) 


< Err 


and E 


e: 


(RE) 


> Et,t- (or reversely), for some given threshold 


We plot in Figure |4ja) (resp. b and c) the error ratio Cr 
( resp. Ce and (^i) obtained for a wide range of Pt (resp. Et 
and Er.r) and several More precisely, the plotted 

stands for the ratio averaged over various RS-to-BS distances 
(Dfe G [600,1000]m), various user rates (TZ G [2,4]bits/ch. use) 
and the two outdoors propagation environments described in 
1 12l Appendix A]. For the purpose of validation, we consider 
a wide cell coverage by fixing the outage requirement Pout to 
0 .1, which is very large for encoded data. 

From Figure HJa), we observe that the error ratio for the 
REA does not exceed 3% when the circuitry consumption is 
not considered = 0). Although the model for E^'^^p^ > 0 

is less accurate, such error increase does not impact at all 
the proposed model for energy consumption, as illustrated in 
Figure Hfb). Indeed, when an error occurs and a user M{x, y) 
is wrongly declared in while it is not (or reversely), we 
have Eb + Er ~ Ed- For the EEA, the error ration does 
not exceed 1.5% and, for the interference approximation, it 
is below 5%, as plotted in Figure Hfc). Furthermore, error are 
mostly located at cell edge and, by considering restricted cell 
coverage (Pout = 0.02), as for the rest of this paper, the error 
ratio for the ICI falls under 2.5%. 

The proposed models for the REA, EEA and ICI thus 
provide very efficient frameworks for performance analysis 
with regards to both accuracy and savings in the simulation 
time. Non-model based simulations were also performed tak¬ 
ing 50 000 samples for the channel gains of each link and 
required several hours for a single relay configuration. By 
comparison, model-based simulations were completed in less 
than 3 seconds. Subsequently, we will refer to ’’simulations” 
for ’’model-based simulations”. 


A. Models validation 

For validation of the proposed model for REA, we account 
for all users M{x,y) located in the simulated (^r,Pt) but 
not declared in ^.4 r, Pt^ , for some given threshold Pt (or re¬ 
versely, M{x, y) is declared inside while it is actually outside). 
This means that, for such user, the effective relaying probabil¬ 
ity, obtained by Monte-Carlo simulations, is s.t. Prtx > Pt but 
the proposed lower-bound gives P^w {x, y) +Pcr {x, y) < Pt 
( or reversely). We define the error ratio ([r as the proportion 
of such erroneous users, i.e. 

^ ff Is^dxdy 
^ ff T^ARdxdy 


B. On the minimal energy consumption per unit area 

As a first step, we do not consider the impact of relays 
in terms of interference and analyze the relay performance 
regarding the EEA only, in the shadowing model used in the 
WINNER project. The Joule-per-bit metric has been widely 
used for energy efficiency analysis. Yet, in practice, a large 
part of the network is primarily providing coverage and does 
not operate at full load, even at peak traffic hours. Due to the 
energy dissipated in circuitry to maintain the network oper¬ 
ational, the energy efficiency can be particularly poor under 
low-traffic loads and restricted coverage Q. To capture the 
aspect of the cell coverage, we consider the maximal energy 



























10 


c 



Px" 

(a) Relay Efficiency Area 



Ex' in mJ 

(b) Energy Efficiency Area 


^ 3 


20 40 60 80 100 

Ex,r in mJ 


(c) Approximation of E ( Ei- 


) 


□ OmJ □ 50mJ □ lOOmJ 


Fig. 4: Validation of the proposed models; error ratio ( for various simulation settings 



consumption Emax that is required to send data at rate TZ to any 
user located within the cell sector, i.e. ,4sector G (-Ae, Emax), 
and divide it by the sector area ,4sector- It is expressed in Joule- 
per-square-meter and denoted as 4/; 

,T, Emax + Eidle . v/3p,2 

4' = --- where Aector = — D^,. (22) 

^sector ^ 

Here, Ey accounts for the transmit energy E as in Eq. 

doll, the RF amplifier coefficients and rjB and the energy 
E(dsp) dissipated at the relay for decoding and re-encoding. 
Eidle refers to the other energy offsets described in Section 
IIII-BI for the considered coding scheme. In Figure we plot 
the minimal feasible energy per unit area 4* as a function of 
the cell radius Df,, considering different values for Nr and 
E('Jsp). To do so, for each considered set (Df,, Nr, we 

find the location for the Nr relays which minimizes 4'. Note 
that proceeding this way would not have been possible in a 
reasonable time without using the proposed models. 

Remark: For comparison purpose. Figure |3 also plots the 
value for 4^ achieved when the relay location maximizes the 


metric F (namely ’’Optimal for F” in the hgure). We will come 
back to this point in the following subsection. 

Result 1. The energy offset consumed for decoding and 

re-encoding at the relay station has severe impact on the cell 
energy efficiency, except for very large cell radius, where the 
overall energy is dominated by the RF transmit consumption. 

As an example, at D;, =700m, increasing from 0 to 

50mW leads to a degradation in 4fmin of 11% with Nr=2 and 
of 22% with Nr=2>. 

Result 2. Up to Nr=3, increasing the number of relays per 
cell allows significant energy reduction. Passed this limit, the 
gain provided by additional relays is not sufficient to compen¬ 
sate for the idle energy dissipated to maintain the 

network operational. 

With Nr=2, the minimal value 4fmin is equal to 0.74e- 
7J/m^, and adding one more RS (making Nr = 3) allows a gain 
of 28%, where 4'niin then reduces to 0.58e-7J/m^. However, 
there is no significant performance gain for increasing the 
number of RS from Nr=3 to Nr=4. 

We recall that we have considered as performance metric 
the maximal energy necessary to transmit data at a given 
rate to any user of the cell sector, as given by the EEA, 
rather than simply the average energy consumption of the 
cell. Simulations show that the relay configurations optimal 
for 4^ and for the average do not match. When optimized 
for the average consumption, the relay configuration results 
in a severe energy increase at cell edge (from 10% to 25%), 
meaning that there exists a major performance gap between 
users of the cell center, with strong channel, and cell-edge 
users, with weak channel. Thus, the average optimization 
criteria does not provide fairness as does our proposed Energy 
Efficiency Area and related metric 4^. 

C. A new energy-interference trade-off on relay deployment 

We now investigate the network performance in terms of 
both energy and interference, by using the F-metric. In Figure 
|6] we plot the maximal feasible Fmax as a function of the cell 
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radius D^, for different values of and To do so, for 

each considered considered set (D^, we find the 

location for the N,. relays which maximizes T. 

First, for a given number of relays, we observe that Fmax 
is poor for small-size cells (D{, <600m) and is significantly 
affected by For example, at D{, <550m, Fmax is divided 
by 2 when E^'^^Pi is increased from 0 to 50mW. On the 
contrary, for wider cells, the overall energy consumption is 
dominated by the RF transmit energy and the interference issue 
is relaxed due to distance. The impact of E^'^^Pi is minor and 
Tmax increases. For the sake of clarity, the case E*^'^^Pi=0mJ 
has been plotted for Nr=2 only, but results do not change for 
larger number of relays. 

Result 3. Except for very short cell radius, the energy offset 
E('^^p), consumed for decoding and re-encoding at the relay 
station, has little impact on the energy-to-interference ratio F. 

Second, as depicted in Figure |6l increasing the number 
of relays generally improves Fmax- We have shown in the 
previous subsection that relay configurations with Nr=3 and 
Nr=4- provide around the same minimal energy consumption 
per unit area T'min- On the contrary, when accounting for 
the interference generated by relays, the case Nr=4- largely 
outperforms Nr=3 and, for example, Fmax is increased from 
2.46 to 3.66 at Dt,=800m. While the transmit energy gain 
achieved by increasing the number of relays from 3 to 4 is 
just enough to compensate for the additional idle energy 
(which is consumed whether or not data is transmitted and 
affects both wcain and 'kmin), it is largely beneficial for the 
neighboring cells (reduced uloss)- 

Result 4. Accounting for the interference generated by relays, 
deploying many relay stations potentially closer to cell edge 
but transmitting at lower power is more efficient than deploy¬ 
ing few relay stations far from cell edge but serving a large 
part of the cell. 

We now compare the performance achieved when the relay 
configuration is optimized either for T' or F and propose 



Fig. 7; Optimal relay positions {Nr = 2, E^^) = 50mW) 


guidelines for efficient relay deployment. To do so. Figure |5] 
also depicts the value for 'F obtained when the relay location 
is optimized for F (namely, ’’Optimal for F” in the figure) 
and reversely. Figure |6] depicts the value for F achieved by a 
location optimized for 4* (’’Optimal for T*”). A first important 
remark is that both relay configurations are essentially distinct 
and provide notably different performance. To illustrate the 
gap between such deployment options, we plot in Figure |7] 
the relay configurations optimal for F and for T' with Nr = 2 
and various D;,. Moreover, we observe from Figure |6] that the 
value for F achieved with an energy-efficient relay deployment 
(optimal for T') is below 1 for 1000m < Df, and Nr > 3, 
meaning that the network performance is actually degraded. 

Result 5. Energy-efficient relay deployment does not nec¬ 
essarily lead to interference reduction and reversely, an 
interference-aware location is suboptimal for the cell energy 
consumption. 

Based on the above results, we propose a guideline for 
efficient relay deployment regarding both T' and F. First, for 
short cell radius (550 < < 700m), deploying two relay 

stations per sector located to minimize the energy consumption 
per unit area T' can be considered as the best option, F remains 
above 1, meaning that the overall network performance is not 
degraded. 

Second, for wider cell size (700m < □&), deploying four 
relay stations per sector provides the optimal results for both F 
and T*. However, current cellular networks are already reach¬ 
ing saturation and negotiating new site agreement for antenna 
deployment is getting ever harder for cellular operators. Thus, 
we argue that considering Nr = 3 may actually be the best 
practical choice. When 700m < D;, < 1000m, optimizing the 
relay deployment for F (resp. T*) does not degrade too much 
the performance in 'F (resp. F), such that both deployment 
options can be considered. However, for 700m < Df,, a 
deployment optimized for 4^ should not be considered since 
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the overall network performance is degraded (F < 1). 


D. Impact of the relay coding scheme 

Up to now, we have shown that a relay deployment can 
be energy-efficient at the scale of a single cell (measured 
by but without necessarily being efficient at a larger 
scale (measured by F). We now investigate the performance 
achieved by the energy-optimized relaying schemes described 
in Section IIII-BI Regarding the F-metric, EO-PDF maximizes 
the energy gain ucain^ while decreasing at the same time 
the energy loss wloss experienced by neighboring cells. As 
detailed in Olhll . the energy transmitted by EO-PDF is more 
uniformly spread over the two transmission phases and over 
both the direct and relaying links, reducing at same time the 
power peaks causing high interference. The IR-PDF scheme 
minimized the use of the relay station by transmitting the 
most data possible via the direct link. The energy loss uloss 
is minimized, but in return, the energy gain ucain is reduced. 
We recall that the circuitry consumption of such partial DF 


schemes is expressed as + Ep'^f^’^'', where 

is an additional offset accounting for their increased 
complexity. 

1) Objective and simulation settings: For analysis, we con¬ 
sider a cell sector aided by two RS only (Nr = 2), with low, 
medium and maximal cell radius (Df, S {600,800,975m}). 
Such configuration provides suboptimal performance in both 
'k and F, compared to a configuration with more RS, but it 
offers valuable infrastructure cost reduction and deployment 
simplicity for a cellular operator. We consider as performance 
basis an energy-efficient relay deployment where both RS are 
located to minimize the energy per unit area 'k consumed by 
two-hop relaying. 

To investigate how optimized relaying schemes can alleviate 
the interference issue, we derive the optimal utilization of 
coding schemes within the cell sector. To do so, we compare 
for each user location the performance achieved by two-hop 
relaying, EO-PDF and IR-PDF and select for each the one 
that increases F. Proceeding this way, we design a spatially- 
optimized utilization of coding schemes and draw a map 
showing the cell areas where to use each coding scheme to 
maximize F. Such map is illustrated in Figure [8] for Db=800m 
and e[,j 7+^ G [0, 50]mJ. 

2) Spatial analysis: With a medium cell radius (Db=800m), 


-(dsp-r) 


-(dsp-r) 


is below 


when the additional circuitry consumption 
lOmJ, EO-PDF outperforms two-hop relaying in almost the 
whole cell area, as illustrated in Figure |8] For higher values of 
EO-PDF does not provide sufficient reduction in the 
transmit energy around the RS and cannot compensate for the 

:50m J 


dissipated energy However, even with Ep^^ 

(i.e. EO-PDF consumed twice as much energy as two-hop 
relaying to process data), EO-PDF still outperforms two-hop 
relaying when the user-to-relay link is weaker. In larger cells 
(Df,=975m), EO-PDF outperforms two-hop relaying for any 
user location and any value of Ep^}’’^^ G [0,50]mJ. The 

additional circuitry consumption Ep^}'’^^ has only marginal 
effect on F. In smaller cells (Db=600m), the overall energy 


-(dsp-r). 



Fig. 8; Optimal utilization of coding schemes in a cell sector 



Fig. 9: Fmax with optimal utilization of coding schemes 


consumption is dominated by the circuitry consumption. Ex¬ 
cept from the case Ep'^}’’^ = the EO-PDF scheme 

improves the cell performance only if used very far from RS. 

Also note that the IR-PDF scheme can only reach the same 
performance as other coding schemes but without outperform¬ 
ing them. The corresponding cell areas are plotted in green. 

3) Coding schemes and relay deployment: In the following, 
we denote ’’combination EO-PDF / 2Hop” as the spatially- 
optimized utilization of coding schemes, as previously de¬ 
scribed and illustrated in Figure [8] and plot in Figure |9] the 
maximal energy-to-interference ratio Fmax achieved by this 
combination and by two-hop relaying only, with Nr = 2 and 
Nr = 3. 

Remark: For some data points, two-hop relaying only 
outperforms the combination EO-PDF / 2Hop. Indeed, such 
combination is used with an energy-efficient deployment (op¬ 
timized for 'k), which is sub-optimal for F. On the contrary, 
the plotted performance of two-hop relaying is obtained with 
an interference-aware deployment (optimized for F) and can 
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be understood as the maximum T reached by two-hop relaying 
in the best feasible relay conhguration. 

Result 6. In a small-size cell, the energy-optimized partial 
decode-forward scheme (EO-PDF) is severely affected by its 
increased circuitry consumption such that the optimal combi¬ 
nation EO-PDF / 2Hop does not provide much performance 
enhancement when >30mJ. 

For example, when Db=600m and =50mJ, the EO- 

PDF / 2Hop combination achieves r=1.33, which outperforms 
two-hop relaying with same relay location (r=1.19) but is 
below the performance of two-hop relaying when used with 
an interference-aware location (r=1.67). 

Result 7. When the additional circuitry consumption Ep^^’’^^ 
is low, or when the cell radius is wide, the combination EO- 
PDF / 2Hop approaches and even outperforms the maximal 
Pmax achieved by two-hop relaying, even with a relay location 
suboptimal for F. It also allows a reduction of the number of 
relays per sector for the same, or even better, ratio F. 

As observed in Figure |9] the maximal Fmax achieved by the 
combination EO-PDF / 2Hop is higher than two-hop relaying 
with Nj. = 3, for almost any cell size and any additional 
circuitry consumption. It even outperforms the case = 4 
for wide coverage extension (Db=975m). 

VIII. Conclusion 

We have highlighted a new trade-off on relay deployment 
for cellular networks that balances system energy efficiency 
and performance loss experienced by neighboring users due to 
the additional interference generated by relays. To this end, we 
first formulated a spatial definition of the relay efficiency and 
proposed three tractable models allowing meaningful analysis 
without requiring time-consuming simulations. Next, we ana¬ 
lyzed the correlative impact of the circuitry consumption, the 
location and number of relays as well as the relaying coding 
scheme on the network performance. By signihcantly reducing 
the transmit power peaks, energy-optimized coding schemes 
alleviate the interference issue, and by performing well even 
with suboptimal relay location or reduced number of relays, 
they offer valuable deployment flexibility. 
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Regarding the probabilities of lines (a) and (c), the condition 
Et + Er < Ed necessarily holds given energy constraints and 
these two probabilities can be readily computed in closed- 
form. Second, the probabilities of lines (b) and (d), denoted 
p(^^ and P('^^ respectively, can only be expressed in integral 
form, but are respectively upper-bounded by 

p(^) < P -f <Edf^ <Esyj E^"’^ < Er^j 

P(‘^) < P < Ed < E^^^ -f 

n E^'"^ <Es + Er< E^^^ -b E^^^ ) 

Plugging these upper-bounds into the expression for 
P {Ef, + Er < Ed), we obtain the lower-bound P[J^ given in 
Eq. ®. Next, we have; 

p(2) = P <Ed< E^'^^ n Eb + Er< 

+ p <Eb + Er<Ed< E^"’^ n 

Note that the energy Er consumed by BS to transmit data 
to RS is generally low thanks to strong channel conditions. 
Thereby, the probability in second line approaches 0 and P^^^ 
can be tightly lower-bounded by Pj^^^ of Eq. (fTOl l. We now 
show that Plow = Plow + approximated by a 

closed-form expression. First, P^w, and thus Pcr and Pcd, 
are computed using the cumulative distribution function $ of 
the standard normal distribution; 

p [Ek < ^ - (24) 


Appendix A 

Proof of Lemma[T1 Lower-bound for Per 
The probability Per is expressed as follows; 


r(m)p(m)p(m) 


J J J ^ (7?b + Er < Ed) dEddEbdEr (23) 


0 0 0 


Among the wide possibilities for lower bounds, we aim at 
removing the power constraints which condition Eb + Er < 
Ed- Hence, we decompose Per into elementary probabilities 
that discard the triple integral. We get Per = P^^^ -b P^^^ with 


p(i) ^p(^Eb + Er<Ed< 

= P <Ed< E^"' n Er< E^"^ n Ei, + Er< Ed^ . 


where /ifc and ak are given in Eq. dSll. Even if $ is written as an 
integral, it is widely available in scientihc tools through well- 
known tables, such that its computation does not imply much 
complexity and can be considered as closed-form. In addition. 
Plow requires the computation of P ^Eb + Er < Eg^^ + 

and ¥{Eb + Er < Ed) = ]? ^ + E^'^ ^ , both 

of which involve the sum of two log-normal random variables. 
Such distributions do not have a closed-form expression, but 


have been extensively explored in the literature 112511261] . In this 


work, we consider the Fenton-Wilkinson approach and approx¬ 
imate these sum distributions by log-normal random variables. 
Eb + Er is approximated by Eb+r ^ lo^Af (ub+r, o'?,...), 
where ^b+r and ab^r computed as given in 
12)]. Similar computation can be performed for E^^ 


: ^'b+rj 

Eq. (9- 
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taking into account the correlation coefficient between 

and Ei°'>^. 

0 Si, ' Sr 

Consequently, using the Fenton-Wilkinson approach, we 
have decomposed Piow(a;, y) into elementary probabilities that 
can be computed in closed-form using Eq. (l24l i. thus avoiding 
the computation of a triple integral for each possible user 
location M{x,y) within the cell. 


Appendix B 

Proof of Lemma|2J Bound for PerE [Ei,+r \ Cer] 

Here, we use the decomposition for Pio* that has been 
proposed in Lemma [T] 


^ErE [i?b+r I Cer] < El + E2 


with El = P^^^^E 
E2 = P^^^E 


Eb+r \ Eb + Er < Ed < E 


(m) 
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and 


Eb+r I E^'"^ <Ed< E^'" n Eb + Er< E;-" 


-(m) 




where g is given by Eq. (Ull. With regards to Ei, we come 
back the integral form. Denoting fk the probability density 
function of Ek, we get; 
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